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ABSTRACT 


The technique of reconstruction tomography used widely 
in diagnostic radiology has been successfully adapted to measure 
density in bubbly air— water flows. The ^convolution back- 
projection' algorithm^ in conjunction with the Ramachandran- 
Lakshminaraynan filter has been used for reconstruction in the 
form of a density distribution of a cross— section. The process of 
data acquisition is governed by photon statistics, which follow 
the laws of Poisson random variables. 

In this study, a Poisson error has been imposed on the 
projection data and the consequent effects on the reconstructed 
profiles has been studied. It has been observed that the process 
of reconstruction is more vulnerable to uncertainities in the 
region very close to the centre of the pipe. Also, the need for a 
gamma— ray source stronger than 13 mCi is reflected in the results. 
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CHAPTER 1 


TMTRODIJCTION 

Compute r i 7Pd tomography, which revolutionized the area 
o-f medical imaging, is being increasingly used in a variety of 
non-medical areas such as nondestructive testing (NDT) and 
two-phase flow measurement. One of the major problems confronting 
a nuclear engineer is the measurement of void fraction during a 
"loss of coolant accident" (LOCA). Tomographic techniques are very 
useful in such instances due to their accuracy and reliability. 

In many appl i cat i ons , in which measurement of 

cross-sectional distribution of any property is required, we can 
make only indirect measurements by probing the object with 
invisible, penetrating radiation, and then interpret these 
measuremen ts . Often, this measured data is not directly 

inter pretable , but is related to the cross-sectional distribulior 
(of the relevant property) in a known way. The general aim of all 
image reconstruction procedures is to process the data to form i 
cross-sectional image, thus facilitating the interpretation of the 
measurements. To determine the density of the object under test, 
various strip integrals (of a parameter like the attenuatioi 
coefficient ) corresponding to a particular angle of view an 
taken. This set of strip integrals is called a projection of th( 
object. Given a number of such projections at different angles o- 
view, the estimation of the corresponding distribution (o 
attenuation coefficient or another parameter ) within the objec 
is the basic problem of image reconstruction from projections 



Computerized Axial Tomography (CAT) is the most significant 
application to date, of this technique. 

The mathematical principles, which form the basis of 
this technique, are attributed to Radon C13, who showed that any 
srbitrer^ -function can be reconstructed provided c.ii of its line 
integrals are known. However, the method of Radon could not be 
implemented due to the inherent mathematical comple;jities. 
Bracewell reported an application of CT in radio astronomy CEU and 
Cormack LEI! derived an inversion formula which was closer to being 
implementable compared to Radon's solution. Bracewell and Riddle 
II2I! and Ramachandran and Lakshtriinarayan C4I] showed how 
computations involved in the CT methodology could be accelerated 
by adopting the use of convolutions. However, it was almost 50 
years after the publication of Radon's paper that a CT machine was 
pioneered by Hounsfield C5I1, in the year 1970. 

The applications of tomography can be found in such 
diverse and wide ranging areas like radio-astronomy CEO and 
electron-microscopy C40. Butchers have used CT scanners to 
estimate the quality of meat prior to slaughter C60 . Use of CT 
scanners to inspect wooden poles used in power transmission has 
also been reported C7D. 

The present study is an effort towards studying the 
effect of the errors in the data collection by a gamma— ray 
counting set-up. The physical laws of nuclear decay are such that 
their intrinsic statistical variations have a Poisson 
distribution. This implies that the projection data used in 
reconstructing the image has an error term inherent to it. In this 
study, this error term has been replicated by incorporating 



pseudo-random numbers in the projection data and the 
effects on the reconstructed profiles investigated. 


consequent 



CHAPTER a 


PRELIMINARIES 

The interaction of radiation with matter is through 
scattering and absorption. Absorption of photons , during their 
passage through matter results in attenuation of the beam. This 
property is Used in reconstruction of the density profiles of the 
two— phase flow. 

Single beam mono— ene rgeti c radiation phenomenon in a 
plane is represented by 

N = N^eKpj^ -S j. (1) 


The value of pi is characteristic of the material and 
also depends on the energy of the incident radiation. Eq (1> 

considers /j to be a two-dimensional function of positionr as the 
path of radiation is assumed to be restricted to a plane. For 

monochromatic radiation the energy dependence of is not of 

concern. For a cross-section of interest, having non-uniform 

distribution of ’pi', equation (1) reduces to 

p = ln(N /N) = f pi{r,<s!>>dl (2) 

o 

c 


Radon Cm showed that it is possible to recover ju from a 
set of several p-values measured along various chords, c. If 

desired, these ju-values can be suitably calibrated to give the 

✓ 

density values- This property is utilized in estimating the void- 
fraction in a two— phase flow. 



2.1 DATA COLLECTION MODES 


Data collection in tomography defines to a considerable 
extent, the speed of the reconstruction and the complexity of the 
reconstruction algorithm. The main modes in data collection are 
the parallel beam-geometry and the fan-beam geometry CSH. 

2.1.1 PARALLEL-BEAM GEOMETRY 


In this mode there are several pairs of radiation source 
and detector systems which scan the object completely. As shown in 
Fig.1, the source - detector pairs are spaced uniformly and the 
object to be imaged is stationed on a rotating table to give 
different values of The line SD represents the path of the data 
ray. The perpendicular distance of the ray from the origin is 
denoted by Several SD pairs collect the data p for a given &. 
This set of p is known as a ’projection*, which are collected for 
different views ( at different D >. This data is denoted by 
p(s ;&) . 

2. 1 . 2 FAN-BEAM GEOMETRY 

In this mode, a single source is viewed by several 
detectors, simultaneously. The source angle is denoted by <y ant 
the detector angle by f3 . The readings are taken at differen 
values of f3 to get the data, denoted by g(<y,/?) or h(X,/^>. Here, 5 
is the perpendicular distance from the origin to the particula 
ray. This configuration, depicted in Fig.S, is widely used i 



S - Source 
D -Detector 



RAY = (s, 0) 
f = f(r,cj)) 


Fig. 1 Parallel beam collection geometry. 




Fig. 2 Fan -beam data collection geometry. 
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tomography for medical purposes. 


a. a MATHEMATICAL FORMULAE FOR TOMCXjRAPHIC INVEKSIQN 

Tomographic inversion methods are generally categorised 
as the ’transform methods’ and the 'algebraic reconstruction 
techniques’. The transform methods are based on mathematical 
formulae while the latter are iterative computational techniques. 
In the following sections, "transform" methods are briefly 
summarised for the two data collection geometries. 

2.2.1 PARALLEL BEAM FORMULATION 

The reconstruction process of a two-dimensional function 
from its projections (line integrals ) involves Fourier transforms 
for parallel beams and Hilbert transforms for divergent beam data. 
A very convenient transformation from FBG to PBG avoids the use of 
the Hilbert transform. The ’ central— si i ce ’ theoreiii (also called 
the projection theorem ) states that the one dimensional Fourier 
transform of the projection data, p(sj&), with respect to the 
first variable s,is equal to the two-dimensional Fourier transform 
of the object function f. Mathematically, 

f(Fjd) = p(F;^) (3) 

Taking the inverse Fourier transform of (3) we get 

n CD 

f(r, 46 )=J J p(F;d)expC-i 2 rTFr cos (d-<ii. ) 3 j F j dFd& (4) 


o 


“00 



For a given & 


R 

p(F;6^)= J p(s;d)expC-iSnFs3ds. (5) 

-R 

Equation (4) requires continuous projection data for all 
values of s and For computational feasibility, a filter 

function is incorporated in the Fourier domain, to enable a finite 
cutoff frequency. This necessity arises because the limits on R 
vary from -co to -Ko, which introduces divergence. This filter 
function, W(F), vanishes for ]Fj greater than A, the cut-off 
frequency. Eq (4) is transformed to 

n CO 

f{r,<ii»)= J J p(Fjd)expC-i2rrFr cos ) 3W (F > [ F | dFd^ (6) 

O -CO 

The reconstruction is approximate as all the higher 
frequencies have been eliminated. Besides, as per the sampling 
theorem, the cut-off frequency and the sampling interval in the 
spatial domain are related by 

A > C1/(2As)3 (7) 


where As is the spacing between the rays. 

The band-limiting filter, introduced in electron 
micrography by Ramachandran & Lakshminarayan C43 is given by 


r ^r\F\<fi 

W(F)= <8 

t 0,|F[>A 


The sine filter used by Shepp & Logan C93 is given by 



W(F) = 


sinC{rrF/HA)0/(nF/2A) 


1F|<A 

B I— i -w 


(9) 



The inversion formula for the FBG case was first derived 
fay Herman and Naparstek (1978> and is given fay, 

^ zn B 

f(r,(ji') = (1/4n^) / i' ( 1 /sin C<> ’-o- ) ) D 9 ( 0 *,/?) ck/dfi , (10) 

O -B 

where, D = distance of source from reference origin, 

<y = angle of the data— ray with the reference ray, 

B = angle of the extreme data rays, 
ft = source position, 

gic^fft} = data for the ray represented fay (o*,/?), 

Dv = d/U) C 3, 

o-' = tan *'C ( r cos (/3-<?^ ) ) / (D+rsin ) ) □ , 

U = Lircos{ft-4,)y^ + (D +rsiT\{ft-^}f 
and the remaining variables are same as for the PBG case (see 
Fig.1). In the above equations U is the distance of the radiatior 
source from (r,<}l-) ,the point being reconstructed ,and £>■’ is thi 
angular displacement of the particular data ray passing througl 
that point (r,<ji). A major difference ( compared to Eq.{6>, whici 
derives the reconstruction formula for parallel beam geometry ) i 
the computation of partial derivatives of the data, g. 

For FBG instead of using the formula derived by Herman 
Naparstek, we used the transformation C10Ii, 



1 


X 


= r^+tan”* (X/D) 

1 + (X/D>^ 

That is, 

h(X,/^)=hf , /9=&-sin~^ { 1/D> 1 s p(l,^) (11) 

y 5 

where p denotes the source angle and X denotes the 
perpendicular distance from the origin to the particular ray from 
the source. In Fig. 3, this conversion has been graphically 
depi cted. 

2. 3 CONVOLUTION BACK-PRQJECTIQN ALCSORITHM 

Reconstruction algorithms based on the inversion 
formulae are called ’transform methods* 181. They are slow and 
require accurate interpolating schemes in computing the 
two-dimensional inverse Fourier transformi. The introduction of 
convolution 12, 41 eliminated the use of the Fourier transform 
and its inversion. 

For PBG , equation (4) can be written as 

n R 

f(r,<j^)= J J p(sjd)q(s'-s) ds d&^ (12! 

o -R 

A 

q(s)= J W(F) jFj exp(i2nFs)dF (13! 

-A 

and s'=rcos(^-^) 

Here q is known as the convolving function and is th 
inverse Fourier transform of the function, W(F>jF|, where W(F) i 
the filter function. The inner integral of eq (S) is a convolutio 




and the outer integral is called the back— pro j e ct ion . The CBP 
algorithm is the most widely used method of reconstruction 
employed by commercial CT scanners used in the area of medical 
imaging . 

2,4 AN OVERVIEW OF PHOTON STATISTICS 


The quality of reconstruction is affected by the errors 
inherent to the data collection process itself. It is further 
affected by the discretization of the problem ( that is, because, 
the number of projections is finite ), and due to the errors by 
numerical computation. In this work, the errors due to photon 
statistics have been studied. 

A very basic limitation to the accuracy of measurements 
taken in CT is the statistical nature of photon ( whether gamma 
ray or x-ray ) production, their interaction with matter and 
photon detection C11I. In gamma ray tomography, the rays are 
generally monochromatic, or at most, two distinct wave lengths are 
present. Due to this, the problems of uneven detector response to 
the detector, and of beam hardening, which plague x-ray tomography 
are mitigated II1S3. 

If we consider that all the photons emitted by the 
source in a unit period of time in the direction of the detector 
are detected, the possible counts ( which are non— negative 
integers) give rise to a discrete random variable denoted by Y. It 
can be further shown that the probability of occurrence of Y, at a 
specific value y is denoted by 

— exp(— A)A^/yi 


(14) 



where A is a fixed real number 


Eq (17) is called the Poisson probability law, and Y 
which satisfies this law is called a Poisson random variable with 
parameter A. The main properties of this variable are as follows i 

a) its mean is A, 

b) its standard deviation is Va 

c) it behaves normally if A is large (>100). 

To estimate A, which is the average number of photons 
emitted per unit time by a stable source (gamma ray or x— ray ) in 
the direction of the detector, all the photons reaching the 
detector will have to be counted. A will be estimated fay the count 
of the number of photons during a particular period of unit time 
(i.e. by a sample of the random variable). Poisson statistics 
imply that if we increase the sample by a factor of N, we reduce 
the size of the li^ error by a factor of VN, where is the 
standard deviation. 

The interaction of photons with matter also can be 
modeled as a discrete random variable. A photon leaving the source 
in the direction of the detector, will reach the detector (without 
being absorbed or scattered) with a fixed probability v^. This 
probability depends on the energy of the photon and the material 
lying on the line between the source and the detector, yf is called 
the transmittance of the material ( along that particular line ), 
at that particular energy. Of the photons which leave the source 
in the direction of the detector, a fraction v' will eventually 
reach the detector, the rest being absorbed or scattered. The 
photons on reaching the detector are counted with an efficiency £ y 



which is called the efficiency of the detector- 

Thus the number of photons which actually reach the 
detector without being absorbed or scattered, and are counted by 
the detector is a sample of the random Poisson variable with 
parameter . 

As discussed earlier, 

- In (N/N )= r ju(r,<^>dl a* m 

Cj j 

c 

where m is the ray sum and is the sample of a random 
variable such that. 


ju + ln(N/N ) < S 

M O 


where 


and is 
of the detector, 

A is the 

OL 

measurement , 

is the 

OL 

£ is the 

<dL 


S {<^ A V-f ) 

a OL OL a 

the fraction of photons leaving in 

number of photons emitted during the 

transmittance of the material and 
efficiency of the detector- 


( 15 ) 
direction 

period of 



CHAPTER 3 


PROGRAM IMPLEMENTATION AMD DESCRIPTION 

There are various programs developed, in the course of 
this thesis work, which simulate the projected data, reconstruct 
images from actual or from simulated data and those which are used 
just to write or feed the data in a specified form. These programs 
have been written in FORTRAN. The program for simulating an error 
governed by Poisson statistics have been written in TURBO-PASCAL, 
by using the random number generation function available. 

3.1 COMPUTER IMPLEMENTATION OF THE CBP ALGORITHM 


The problem of reconstruction from projections is stated 
as follows s 


Given p(s;P),find f(x,y>. 

In practice, the problem may be stated as s 
Given discrete projection data in the form of estimates 
of p for a finite number of rays, find a E— D distribution, which 
is a reconstructed estimate of the unknown object. 

In the case when p is sampled uniformly in both s and P, 
for N angles Ad apart, with each view having M equispaced rays As 
apart, we define C83 


M^=(M-1)/E 1 
M =-(M-1)/SJ 


M odd 


(M/E)-11 
M"= - M/S J 


even 


In order to ensure that the collection of rays specified 


S5 



by 

CCiriAs^nA^) : M~< 1:£n:S:N>> 

covers the unit circley we have, 

A^=rT/N and As =1/M^ 

A reconstruct ion algorithm which can be implemented on 
the digital computer is required to evaluate f^(kAx,lAy)y which is 
a band limited approximation of the function to be reconstructed. 
Here K <k:^K^ and where k and 1 are the positions of the 

co-ordinates of the image pixel. The definition of their upper and 
lower limits is similar to that of Thus the projected data 

from N views and M rays is to be used to construct an image of KxL 
pixels. In this particular case ,the image is composed of 21x21 
pixels frorri 18 views each having 21 rays. The back— projection 
integral is evaluated as follows : 

w 

f ( kAx y lAy )ii£A0‘ T]p(I^AiLXCOS0+lAysin^ (16) 

B n n ri 

For each angle , the convolved values of p(s",^ ) for 

ri ri 

the KxL values of s' We can either have a separate convolution for 
every S' with the actual value of q(s-'-iTA.s) at that point or we 
can evaluate ’p(mAs,d ) only within the specified limits of m and 

Ft 

then use interpolation. The latter approach is much faster anc 

cheaper. These operations are represented by 

• ^ 

M 

p (nrAsy^ )i!£As !r p(nAsy^ ) q((m— miAs)y (17 

YH Ami r ^ 

m = M 

m < 

p {s',^ )s£A^ P (mAs^ ) 1(3— md s) 

X 1"^ a Fi 


( 18 ) 



where Ks) is an interpolating function 


A 


1 inear 


interpolating function, say I (s), corresponding 

la 

interpolation between adjacent samples is 


I (s) = 


i Or 


! s!<As 
j s I >As 


to linear 


Using the above formulae , a program was written in 
FORTRAN to implement the CBP algorithm. Another program was 
written to simulate data for some standard functions Clike f<r)= 
1, r, exp(r), exp(— r), exp(3r) etc where r is the distance from 
the origin!. 

The data generated by the latter program was used to 
check the efficacy of the CBP program. Certain modifications were 
made in the CBP program to accommodate data obtained in the FBG. 
The mathematical basis of these changes is the FBG to PBG 
transformation mentioned earlier. 

A calibration curve Cdensity <p> versus attenuation- 
coefficient was plotted using reconstruction results for 
water, pine, walnut and air (whose densities are known), after 
necessary corrections for the plexi— glass piping had been made. 
The CTN (Computer tomography numbers) for various two— phase flow 
situations were obtained. With appropriate normalization, the CTN 
and }-t is identical. This normalization procedure, too, was 
included in the program. With the help of the calibration curve an 
idea about the point wise distribution of ct is obtained. The 
density p and the void-fraction ot, of a flow, are related fay the 
expression 


£S= 1-p 


assuming the density of steam to be negligible- 



To estimate the sensitivity of the CBP algorithm to the 
errors inherent to the data collected^ an error with a Poisson 
distribution was deliberately introduced into the projection- 
data« This was done by using the random number generation function 
available in TURBO-PASCAL. 

If N was the original count ( pro j e ct i on— data ) ^ an error, 
which was a random number between — /~T 4 and was added to the 

original data- 

fs} = }h 4 4- random number 

poi. 

where N is the modified data 

pOL 

Errors -For and 3ty deviation (where <y is standard 
deviation ) were also generated, by simply multiplying the random 
number generated by a factor of S.O and 3.0 respectively. 

3. Z PROGRAM FOR DATA SIMULATION 

This program (see Appendix- B) computes the 
line-integrals of a given function between specified limits (in a 
bounded region). This program computes the line— integrals along 
parallel lines at equal distances from each other, for a number of 
views. The various steps of the program are as follows :- 

1) The parameters of the program like the upper limit 
(B), the lower limit (A), the number of divisions (D), the slope 
(M) of the line and its intercept (C) are read from a data file. 

S) The given function is incorporated in an appropriate 
form into the program. 



3) The line integral of the above function is computed 
by using Simpson's one— third rule. Since it is path dependent, the 
integration rule has to be appropriately modified. A subroutine 
has to be incor porated, to take care of the condition when the 
line IS parallel to the Y— a:<is and its slope is infinite. 

4) This whole process is repeated as many times as the 
number of views desired. The value of the intercept ’C’ is changed 
during each iteration to ensure that the parallel lines along 
which the integrals are computed are equidistant from each other. 

This program can be modified to simulate the data for a 
fan-beam geometry also. 

3. 3 PROGRAM FOR RECONSTRUCT! MG THE IMAGE FROM DATA 

This program (see Appendix C) can reconstruct an image 
from either simulated data of Section 3.E, or experimental data. A 
program was written to accept simulated data for the parallel beam 
georrietry and use it to reconstruct the image. Images of several 
standard functions were reconstructed within acceptable error 
1 imi ts . 

This program was modified to accept data from the 
fan— beam configuration, by incorporating the transformation 
mentioned by Kwoh et al C10Z1. The salient feature of this prograrrt: 
not found in C9II (which has reconstruction program for PBG ), was 
to display the reconstructed image as every view was convolved and 
back projected. 

A suitable image quality index would be helpful in 
quantitatively estimating the accuracy of image reconstruction. 



reconstruction 


Two such indices are 1 and 1 C133. The process of 
can be terminated if the error of reconstruction reaches a certain 
minimum value - Another index which has shown promising results is 
the " fractal dimension *' C14I]. The actual program details are as 
f q 1 lows : 

1) The parameters of the program like NDIV. (the number 
of rays in a view ) and NVIEW (the nurriber of views ) are read in 
from a data file. For fan— beam data the corresponding parameters 
will be the source to detector distance (D), the radius of the 
object (RAD), the number of projections in a scan (IMRAY) and the 
number of scans (NVIEW). 

2) The projected data are read in from a data file. If 
the data values are for a fan— beam geometry then the 
transf ortTiation of co-ordinates from fan-beam geometry to a 
parallel beami geometry was done. 

This transformation enables us to use the CBP algorithm, 
for fan-beam data, without resorting to the Hilbert transform. 

3) The Ramachandran-Lakshminarayan filter was used in 
this program. The numerical values had been computed and stored in 
a file. These values were merely read from the file into an 
array. It is at this stage that the convolution is carried out. 

Convolution is especially easy if the values to be 
convolved are stored in an array. 

It is often found that, during back-projection of data, 
the convolved values do not exactly pass through the point being 
reconstructed. Due to this interpolation has to be used. Since 
linear interpolation is fast and it gives good results, it has 
been used in this particular implementation. 



4)The super imposition ( backpro j e ct i on ) of these 
interpolated values is carried out. Some conveniences like display 
of point values as they are reconstructed have been added to the 
program. These reconstructed values are the LITF values. 

A sample output of this program has been displayed in 
Appendix E- 

3. 4 PROGRAM TO INCORPORATE RANDOM ERRORS INTO THE DATA 

This program (see Appendix D) has been written in 
TURBO— PASCAL. It utilises the random number generating function 
available in the TURBO— PASCAL library. 

For a Poisson distribution with mean ’N’, the standard 
deviation is given by ±VH. A random number with a maximum absolute 
value of -i-yN is generated. This quantity is superimposed on the 
experimental data, in effect simulating a Poisson distribution for 
the inherent statistical variations. 




CHAPTER 4. 


RESULTS AND DISCUSSIONS 

In this chapter, the results obtained with the actual 
e:< pe r imental data have been discussed. 

4. 1 DATA USED 

The data used was taken from the study conducted by 
liunshi C15I1. The details are briefly summarized as followsi — 

a) The source of nuclear radiation was 13. S6 mCi of 
Cs— 137. The collimator fixed at the other end was rotated through 

A Ck 

an angle of 50 ( 25 on either side ). 

b) The detector was connected to this, the subsequent 
stages being a preamplifier, an amplifier, a baseline restorer, a 
single channel analyser and a timer. 

c) The data collected for calibration was for water, 
air, walnut and pine. The air— water bubble column was used to 
measure densities fay reconstruction tomography. 

These calibration quantities have been listed in Table 

1 . 

4. 2 DISCUSSION OF THE RESULTS 


The results of the investigations are presented in 
Fig. 4-15 and Tables 3-5. 

The reconstruction results given in Fig. 5-15 have been 
labelled using the following convention 
1) T denotes two-phase flow. 



E) The second numeral (10 


to 


40) 


denotes 


the 


c ross— se c t i onal iy averaged vo id— f ract i on in that particular case. 

3) The alphabet has been used for distinguishing apart 
the different cases, and is connected to the different case 
numbers. No alphabet has been used to denote the first case. 

4) The last numeral (0, 1, E or 3) denotes the magnitude 
of the error, as a multiple of o'. 

Expression (15) in section E.4 shows that the error in 
measurements depends on various different factors and it may be 
reduced by 

1) Increasing the strength of the source or the time 
period of measurement and thus <p^. 

E) Improving the efficiency of the detector. 

In Fig. 4, the relationship between LITF and p (density) 
for known materials has been graphically depicted. It is a linear 
relationship and has been used as a calibration curve for cases 
where LITF is known and p(density) or a(vQid fraction) is unknown. 
Thus, with the help of Fig. 4, further investigations have been 
carried out, regarding the deviation in density or void-fraction, 
by relating it to the change in the value of LITF. 

The cases of two-phase flow investigated are listed in 

Table E. 

The projections for the calibration quantities are 
assumed to be relatively less affected by statistical errors since 
averaging over a large number of readings has been carried out. In 
Fig, 5, the reconstructed profiles For these materials, along the 
largest chord (i.e. the diameter) have been plotted. Fig. 6-9 show 
the reconstructed profiles for case 1 with errors from O to So*. 



These have been reconstructed for <a>=0. 1 to 0.4. For an error 
larger than Icy, the error in reconstruction is very large (up to 
60.0 But, Prom the Poisson statistical curve, we Isnow that 
about 99"/ of the random variables lie within a range of ± 3o- from 
the mean for that particular distribution. Thus, if we have a 
reconstructed profile, from data within a broad band of 'true 
counts t 30 " the reconstructed profile also occupies a 'band'. 
The degree of confidence in the reconstruction band increases as 
we increase the range in which the detector counts may lie. 

Various other observations can be made on the basis of 
Table 3 and the radial profiles plotted in Fig 10-15. It is seen 
that the deviation due to the error incorporated is minimum at the 
edge of tiie reconstruction and increases towards the centre. The 
deviation is maximum at the centre of the reconstruction for all 
the cases. These values of maximum deviation have been listed in 
Tables 3- 5. Table 3 depicts the deviation at the centre as a 
percentage of the central pixel value in the original 
reconstruction (i.e. without any error incorporated ). Table 4 
depicts the differences between original pixel value and the 
values M^ith the error incorporated. Table 5 depicts the deviation 
at the central pixel, for Ea and Ea errors. 

If we consider only the magnitude of the deviation, in 
case 1, <c<>=0. 1 has a deviation of 17.647., which increases to 
21.947. for <cx>=O.H . For <a>=0.3, it reduces to 13.977 and again 
it increases to 34.737 for <c«>=0-4 . 

For case 2, <cit>=0.1 has a deviation of 18.857, which 
reduces to 18.267 for <a>=0.2 -It further reduces to 13.857 foi 
<a>=0.3 and then it increases to an all time high of 42.937 for 



TABLE 1 


CALIBRATION TABLE 


Sr - no 

Mat e r ial 

LITE 

Dens i ty 

1 

AIR 

0.0000 

9 

0-0 g/cm 

2 

\ 

PINE 

0.0833 

0.41 g/cm^ 

3 

WALNUT 

0.1480 

0.73 g/cm® 

4 

WATER 

0.2205 

1.00 g/cm® ! 


TABLE 2 

CASES INVESTIGATED 


Void fraction <o(> 

9 

Density g/cm 

107. 

0.9 

207. 

0.8 

307. 

0.7 

407. 

0.6 











TABLE 3 


TABLE OF DEVIATION AS A PERCENTAGE 


S r - n o 

D E 

V I A T I 0 N 

c%> 

1 

J'C'i>=0. 1 

<oi>=0.2 

<oi>=0.3 

<oi>=0.4 

1 

17.6^ 

-21.94 

13.97 

34.73 

H 

18.85 

-18.26 

13.85 

42.93 

3 

18.98 

-13.96 

-24 . 05 

33.48 

4 

18.21 

-19.20 

-25.77 

1 .73 

5 

-19. 16 

18.06 

-25.95 

31.64 

6 

18.82 

-19.99 

16.84 

41.96 


TABLE 4, 

ABSOLUTE DEVIATION OF DENSITY 


Sr.no 

DEVIATION 

<ot>'=0. 1 

<a>=0.2 

__ 

<« >=0 - 3 

* 

o 

II 

A 

V 


0. 1436 

-0. 1734 

0.0897 

0. 1870 


0.1649 

-0.1374 

0.0854 

0.2059 

3 

0.1595 

-0.0923 

-0. 1402 

0. 1305 


0.1711 

-0. 1422 

-0.1488 

0.0085 


-0.1391 

0.1 128 

-0.1430 

0. 1298 

6 

0.1704 

-0.1474 

0.1131 

0 .-2086 











TABLE 5 


ABSOUn'E DEVIATION OF DENSITY FOR Ho- and So- 
ERRflR FOR CAGE 1 


> 

DEVIATION 

Hey 

30 ' 

0.1 

0.0690 

0.1063 

O.H I 

-0.0696 

-0.1005 

0.3 

0 « 0454 

0.719 

0,4 

0.0868 

0.1323 






ot=0«4 « 

In case 3, «=0.1 has a deviation of 18.987. while ot=O.E 
has a deviation of 13.96X . This increases to 24.057 for oi=0.3 and 
then it increases still further to 33.487. for oi=0.4 . 

In case 4, oi=0.1 has a deviation of 18.217. while ot=0.2 
has a deviation of 19.207. . This increases to 25.777 for ot=0.3 and 
then it reduces to 1.737 for ot=0.4 . 

In case 5, o(=0.1 has a deviation of 19.167. while oi=0.2 
has a deviation of 18.067 . This increases to 25.957 for ot=0.3 and 
then it increases still further to 31.647 for ot=0.4 . 

In case 6, o(=0.1 has a deviation of 18.827 which 
increases to 19-997 for oi=0.2 . This reduces to 16.847 for ot=0.3 
and then it increases to 41.967 for oi=0.4 . 

LITF, which is proportional to the point density, 
decreases logarithmically with *N*, as per eq (2), where ’N’ is 
the number of counts, and is proportional to the number of photons 
reaching the detector. Therefore, as a increases, LITF should 
decrease in magnitude. This is because LITF is proportional to the 
density, it being a measure of , the attenuation coefficient. 
This implies that the magnitude of the deviation should increase 
with increasing a. 

If we consider the magnitude of the deviation for a=0.1 
and 01=0.4, this trend is followed in 5 cases out of 6, case 4 
being the sole exception. If, however, all the values of ot are 
considered, then it is apparent that any explicit relation between 
the magnitude of the deviation and LITF cannot be readily 
inf erred-This may possibly be due to a number of reasons- 


Some of these are • 
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Fig. 6 ; Reconstructed profiles for case 
with 'alpha' = 0.1 and error upto 3 'sigma 








Radial position (cm) 

Fig. 8 : Reconstructed profiles for case 
with 'alpha' = 0.3 and error upto 3 'sigrr 




jj_n p9pnj].suoooy 


44 


Radial position (cm) 

Fig. 9 : Reconstructed profiles for case 1 
with. 'alpha' = 0.4 and error upto 3 'sigma'. 





46 



Fig. 12; Reconstructed profiles for case 















I’Tistead of the line integrals, we have strip integrals 
as the projection data. This is due to the finite source 
dimensions, leading to a finite, non zero beam width. 

Values of LITF are averaged over the cross-section and 
calibrated using known density or void— fraction values, for 
corresponding values of <o(>. 

The detector for photon detection had been collimated to 
accept only those photons which come straight from the source. 
However, perfect collimation is difficult to achieve. This allows 
some scattered radiation to enter the collimator window, 
artificially boosting the counts. Since only straight attenuation 
photons are related to the chordal density, scattered photons 
produce an overestimate of the voi d— f ract ion , i.e an underestimate 
of the density (LITF). These effects due to photon scattering have 
not been taken into account. The increased number of counts due to 
scatter will not result in an overestimate of <oi> if the effect of 
scattered photons is the same along all chordal positions. This is 
a fairly strong assumption and may not be always valid. It has 
possibly resulted in underestimates of density (LITF), partially 
compensating for the greater effects of counting statistics at 
lov/er counts C15I]. 

3) Due to the formation of bubbles, and their consequent 
collapse or liberation, the flow is said to be unsteady in the 
local sense, while being steady as a system since no mass 
accumulation occurs in the system. The bubbles generated are of 
different diameters and have different velocities. Due to these 
processes a "dynamic bias" error is inherent to the system. It is 
analogous to the "patient motion" errors encountered in medical 



imaging and is an experimental error. Its man i f as tat i on is in 
form of streaking artifacts across the tomogram. 


/ 



CHAPTER S 


COMCLUSTOMS AMD RECOMMEMDATIOMS 

II has hnen observed in the course of study that the 
1 inear relationship between LITF and p (density) can be extended 
to got roean^ngful results for measurement of vo i d— f ract i on in a 
two-'-phase flow., [t, is also possible to have LITF to be identical 
lo /Jr th^ a I j on constant of the materia I at that energy , by 
suilal'*1y no f ma I 1 7 ? ng the input data and ensuring that the same 
sysf-of!, of units is ronsi stent during computation- 

riio ei^rnj' in reconstruction, due to the error 
incorporated in the projection data folloivs the predicted trend 
only for certain readings in a given set. This implies that, in 
addition to tire effects of photon statistics, other errors like 
dynaniic bias and counting of scattered photons are inherent to the 
p r o j G c I i o n da t a „ 

Tiia siqn of the maximum deviation is negative in about 
337. of the reconstructions. But the reconstructed profile has no 
preference for ^ign for the individual pixels within a 
r e con 5 t rue 1 1 o!i « It is concluded that the statistical error has no 
proference for sign. It unde res t imates or overestimates the 
reconstructed values with almost the same frequency- 

lo reduce the error in reconstruction due to the 
statistical nature of photon behavior, the practical recourses are 
as Follows s 

1) .incre->se the number of detector counts (photons) by 



II 5 i n ’.f s « ) u r r e of ii i q 1 1 e r s t r e n g I h «i In medical a p p I i c a 1 1 o !i s r t h i «= 
won Iff a hiqbor ladiation dosage , rendering this recourse 

una c cc^ fHrihl e « 

2) Use a 1 OV-* energy source for a longer interval of timey 

eo L lunl the effective** mean counts are recorded due to the effects 
of loiiiporal averaqjng^ This would be acceptable only if the 
dura I ion of fijpasuj r omen t were subs tan t lal ly smaller than the 
duralion of Ihe Fastest transjent in the two- phase flow system™ 

JL IS n I' en t from tFie above r e comiiienda t i ons thst 
further investigations are needed along these lines* E^esideSy the 
study was conducted using only the Ramacljandran— Lakshminar ayan 
filter. In VP s t i qa 1 1 oris for otfier w’lndow-f uncti ons r which are 
possibly Jess •’Mjinerable to statistical s^ariations can be carried 
oul« ffip t omng f apl"i i c inversion has been carried out using limited 
<vi{Oiml of flala„ A larger number of vie’wSy as well as rays in e 
VI ^,,(1 I 1 {' r, 1 j { I, i XI h V? tier results* 

It IS rerommpnded to Follow the courses mentioned above 
i,f) nv/pctiqate the effects of the s tali sti cal errors , inherent to 
lifp dula and c f3n5-<=v:|uen t ly y to reduce these effects. 



REFERENCES 


1 J. Radon, "Uber die Bestimmung von Funktionen durch ihre 
Integ ralwe rte l^ngs gswisser, Maanig-faltigkeiten" , Berichte 
Sichsische Akademie der Wissenschaf ten Leipzig, Math.— Phys. 
Kl., ^ (1917), pp. 262-267. 

2 R . N .Bracewel 1 , A. C. Riddle, "Inversion of Fan Beam Scans in 
Radio Astronomy", Astrophy . J. 15Q , (1967) pp. 427-434. 

3 A-M.Cormack, "Representation of a Function by its Line 
Integrals with some Radiological Applications" J . Appl . Phys . , 
35 (1964), pp. 195-207. 

4 G.N.Ramachandran, A.V.Lakshminarayan, "3-D Reconstruction 
from Radiographs and Electron Micrographs; Application of 
Convolution instead of Fourier Transforms", Proceeding . 
National - Academy . Science , USA, 6S, 1971, pp 2236 -2240. 

5 G.N. Hounsfield, "Computerized Transverse Axial Scanning 
Tomography : Part I, Description of the System" B_r J 
Radiol . , 46 (1973), pp. 1016-1022. 

6 "CT Scans for Butchers", Science T o day , 22 , 3 (198S), p. 26. 

7 W.H. Miller, "The Design of a Portable CAT Scanner for Utility 
Pole Inspection", T rans . Amer i can . Nuclear . Society (1986), 
52 , pp. 350-351. 

8 R .M.Lewitt, "Reconstruction Algorithms: Transform Methods", 
Proc . IEEE , 71, No. 3 (1983), pp. 390-408. 

9 L.A.Shepp and B.F. Logan, "The Fourier Reconstruction of Head 
Section", IEEE Trans . Nucl . Sci . , NS-21 (1974), pp 21-43. 

10 Y.S.Kwoh, I.S.Reed, T.K. Truong and C.M. Chang," 3-D 

reconstruction for diverging X-ray beams ", IEEE Trans . 
Nucl.Sci. , NS-25,No.3, June 1978, pp 1006-1014. 


11 G.T. Herman, "Image Reconstruction from Projections : The 
Fundamentals of Computerized Tomography", Academic Press, New 
York M980) . 

12 A. Mackovski, "Physical Problems of Computerized Tomography", 

P ro ceedinqs of the IEEE , 71 , No . 3 , (1983), pp 373—378. 


P.Munshi, "Error Estimates for the Convolution Backproj ection 
Algorithm in Computerized Tomography", Doctor of Philos 
Thesis, IIT Kanpur, (1989). 


13 



4 V.Bhattr P.Munshij and J . K .Bhattachar j ee , "Analysis of the 

Performance of Medical CT Scanners Using Fractal Theory", 
Bulletin of Radiat ion Protection , 13 , No .E , (1990), pp 3S-40. 

5 P.Munshi, "Two— phase Flow Studies in the Bubbly Flow Regime 
using a Scanning Gamma-Ray Densitometer", Master of Science 
Thesis, Ohio State University, (1979). 



Scan No . 


APPENDIX A 


to 

W 

W 

Pt^ 

CJ 

tp 

Q 



<N) 

CM 

to 


cn 

to 

to 

tO 

r-. 

o 

nt 

to 

xr 

ro 

rH 

r—t 

ro 

XT 

rM 


CO 

CO 

CO 

CO 

CO 

CO 

oo 

CO 

CO 


r'J 

r4 

r i 

IM 

Ol 

CNJ 

CM 

oi 

rj 

in 

■o* 

O 

in 


Rl 

o 

O 

to 

XT 


CO 

to 

lO 


mSM 

o 

to 


to 

r-^ 

oi 

CO 

CO 


Rl 

CO 

CO 

CO 

CO 

<N1 

r-j 

CM 

CM 


lai 

r J 

r j 

CM 

r-J 


Kl 

to 

t''- 

to 

o 

r i 

rO 

to 

o 

to 


O 

to 

XT 

ro 

t D 

MD 

o 

rH 

0-1 

t'- 

t-. 

to 

tO 

to 

O 

O 

r - 

r- 


fi*-^ 

rH 

rH 

rH 

rH 


r~t 

rH 

r-H 

to 

Ot 

rH 

r^ 

O 

IRI 

O 

cn 

cn 

tO 


•ef 


cn 

to 

Rl 

r-- 

to 

rH 

to 


rH 

in 

xr 

xr 

Eh 

lO 

to 

o 

CO 

rsj 

OJ 


rH 

rH 

Bl 

— 

rH 

I — t 

< — r 


CO 

Oi 

CM 

'O- 

■SI 

CO 

ro 

CO 

PM 

O 

O 

CO 

CM 

rH 

Rl 

r-H 

to 

to 

CO 

(M 

to 



ro 

Rl 

CM 

r-H 

XT 

r-^ 


CM 

1 — 1 

rH 

rH 

ESI 

BnI 

r-H 

r-H 

r-H 

to 

CO 

•vt 

HT 

cn 

Rl 

BE 



O 

t 

Ch 

tO 

CO 

^o 

RH 




to 


rO 

rj- 

ro 

CM 

mm 

H 5 H 



r-- 

rH 

CM 

rH 

1 — f 

rH 

ESI 




rH 


cn 

CM 


cn 

Bl 

BE 



to 

to 

•o- 

to 

to 

o 

Bl 

Bl 



f-H 

r— f 

KT 

•M- 

ro 

CM 

Bl 

El 



IM 


CM 

rH 

rH 

rH 

Bl 

EE 



r-H 

cn 

rH 

fO 

CM 

to 

n 

i-H 

Hi 



• 

r- 

CM 

xj- 


EE 

cn 




csa 

xr 

xt 

ro 

(H 

Rl 

cn 




rH 

CM 

rH 

rH 

rH 

Bl 


EE 




(M 

to 

!>. 

rH 

Bl 

CM 

to 


lO 

o 

cn 

to 

to 

n- 

XT 

to 

CM 

O 


rH 

M- 

•O' 

lO 

rH 

O 

o 

CO 

rH 

O 


CM 

rH 

rH 

rH 

rH 



I 


uo 

CO 

to 

in 

CO 


Rl 

o 

rM 

rH 

rH 

to 


LO 



cn 

to 

PM 


to 


tn 

i-H 



c^ 

rH 

O 


CM 

rH 

r-H 

r-H 




IEmI 

t-H 


MD 

to 

f-H 

rH 


MMI 


BE 

O 


CM 

to 

to 

xr 




Bl 

r~4 


to 


ro 

1 — f 




El 

o 


CM 

rH 

rH 



un 


El 



C^J 

cn 

CO 

cn 


WSM 

BE 

oo 

xr 

to 

(M 

to 

xr 

ro 



BE 

O 

cn 

fSj 

to 

XT 

ro 

rH 



BE 

rH 

to I 

04 

CM 

rH 

rH 

rH 



EE 

rH 

r-H 


r^ 

o 

CM 

to 

cn 

xr 

o 

o 

to 


CNJ 


xr 

ro 

cn 

cn 

XT 

CM 

cn 

o 

to 


ro 

rH 

cn 

CO 


rH 

to 


CM 

rH 

rH 

rH 




it-H 

r-i 

to 

CNJ 

XT 

CO 

tO 

rH 

o 

03 

rH 

cn 


to 

-o 

xr 

CM 

o 


XT 

CM 

cn 

CN3 

to 

XT 

to 

rH 

o 

CO 

txr 

rH 

to 

» 

CM 

tH 

rH 

rH 

rH 



r-H 

r-H 


rf 

to 

cn 

to 

CM 

xS- 

CO 

r-H 

r-H 

to 

fO 

vO 

XT 

ro 

O 

CO 

to 

XT 

o 


to 

-tf 

to 

rH 

o 

oo 

rx. 

r-H 

o 


C'J 

fH 

rH 

rH 

tH 



hH 

,-H 

to 

on 

O 

to 


O 

xr 

rH 

cn 

o 


CNi 

to 

XT 


rH 

cn 

CO 

to 

r-H 


to 

XT 

ro 


O 

CO 

rx. 

f— H 

tO 

1 

CM 

rH 

rH 


rH 



r-4 

rH 


\o 

ro 

(M 


CO 

to 

1 — ! 

O 

cn 

o 

rH 


ro 


rH 

r— ^ 

rH 

o 

fO 

rH 

to 

■CT 

ro 


O 

cn 

CO 

CM 

o 

' 

C J 

r-t 

rH 


rH 



r-H 


to 

O 

XT 

cn 

O 

to 

oo 

CO 

to 

CO 


cn 

to 

to 

\o 

to 

to 

to 

o 

pH 

C^l 

-rf 

H- 

to 

rH 

o 

cn 

CO 

rM 

tO 

rH 

C^J 

rH 

rH 

rH 

rH 



r— < 



in 

to 

»H 

tO 

cn 

to 

ro 

to 

to 

to 

oo 

TO 

to 


fM 

o 

rH 

rO 

o 

rH 

o- 

XT 

to 

rH 

o 

o 

cn 

to 

pH 

’ 

CM 

rH 

rH 

rH 

tH 

— t 


rH 

rH 

to 

r- 

CQ 

O 

CO 

ro 

rM 

to 

CM 

XT 

’ 

•cr 

xt 

OO 

r-H 

to 

CO 

cn 

cn 

to 


tzf 

XT 

ro 

r-J 

rH 

o 

cn 

to 

pH, 

rH 

J 

c^i 

rH 

rH 

r-H 

1— H 

r-H 


pH 

rH 


"XT 

ro 

CM 

00 

CO 

ro 

to 

cn 

O 

o 

CD 

CO 

O 

r-- 

rH 


o 

XT 

o 

CM 


XT 

xr 

CM 

OJ 

rH 

r-H 

XT 

CO 

1 

CM 

rH 

1 — 1 

rH 

*-H 

rH 

r-H 

r-H 

rH 

lO 

to 

CO 

cn 


r-H 

CO 

tO 

to 

ro 

' 

rH 

r-H 

tO 

r>. 

tn 

O 

o 

o 

O 

C-1 

ro 

L/1 

xf 

ro 

tn 

ro 

fN4 

to 

03 

r^i 

CM 

rH 

rH 

irH 

r-H 

r— t 

< — t 

r-H 

r-H 


■rt- 

to 

O'! 

t'- 

o 

to 

to 

r-r 

ro 

to 

CO 


lO 

rH 

o 

O 

cn 

XT 

r- 

CM 

O 

to 

to 

to 

lO 

to 

XT 

r'N 

pH 

1 

CM 

rH 

r“C 

rH 


f—i 

r-H 


rH 

lO 

to 

cn 

o 

rH 

XT 

to 

CM 

CM 

o 


vO 

o 

o 

vO 


n 

tO 

CO 

o 

r>> 

CD 

pH 

rH 

O 

o 

o 

o 

o 

o 

Csl 

I 

CNJ 

CM 

CM 

CJ 

CM 

INJ 

oi 

(M 

PM 


‘ Ot 

to 

vO 

r-^ 

xr 

cn 

to 

to 

OO 

O 

tO 


CM 

to 

to 

f" 

CO 

to 

to 

to 

c:*) 

to 

! CO 

CO 

cn 

co 

CO 

CO 

CO 

1 

CM 

C-'J 

f'j 

CM 

C'l 

CM 

CM 

04 

r 4 









Jm' 










to 

r-H « 









t". 

>-H 

5 


Q 

Q 

Q 

Q 

O 


» 

XT 

r 1 





>— • 


CxC 

» 



Ol 

P 

O 

>• 

P 

O j 

o 


cx 

It 

v5 


< 

tfO 





1-5 

tz 




to 

o 

o 

o 

o 




. P 


XT 

xr 

to 

rM 

r-H 



y 










pH 









rd 




















Scan No. 




CO 

PU 

UQ 

c:^ 

ID 

W 

Q 



CO 

^o 

CTi 

-O' 

to 


o 

t « 

oo 

f 

r >. 

CO 


r'O 

CO 

CO 

CO 

OO 

00 



r 1 

P'l 

C'J 

fxi 

(XJ 


LO 


CO 

tNj 

(XJ 

to 




r-. 

O 

•TO 

CO 




00 

CO 

ijO 

CO 


r^i 



04 

rxi 

fxj 




CTv 

t-- 

VO 

■rr 


ir) 


H- 

to 

to 

rj 


rvj 



rv. 

IX. 

t-> 




rH 

r-4 

rH 

t—i 


vn 


to 

O 

to 

lO 




ro 

r-' 

to 

Kt 


r-j 


to 

XT 

xr 

XT 


r-j 


rH 

rH 

rH 





to 

CTl 

o 



O 


to 

t/4 


to 


r^i 


•^T 

to 

t-O 

(XJ 




rH 

1 — 1 

i-H 

rH 


LO 


to 

04 

vO 

o 


* 


o 

<XJ 

to 

VO 


r-- 


xr 

to 

(xa 

rH 


rH 


iH 

rH 

rH 

rH 





(XI 

vO 

to 


to 


CO 

rx 

IX 

cn 


rH 


to 

<xj 

f — ( 

o 




rH 

•-H 


rH 


to 


o 

MD 

to 

o 


* 


CO 

XT 

(XI 

(X4 


CNi 


to 

(NJ 

1 — 1 

O 


rH 


rH 

rH 

rH 

-H 




vO 

to 

o 

to 


o 


CO 

r 4 

C7t 

rx 


rH 


to 

r 4 

O 

cn 




rH 

rH 

iH 



to 



to 

rx 

cn 




r'>- 

o 

tx 

VO 




to 

CNl 

o 

cn 




rH 

rH 

rH 





cn 

vD 

r-. 

(XJ 




CO 

cr^ 

XT 

to 


ji 


to 

r—i 

O 

cn 




1 — ( 


rH 



LT) 


r^j 

^ 

to 

CO 

cn 




cn 

o\ 

XT 

rH 


fVJ 


to 

rH 

O 

cn 




r-i 

I — 1 

r-* 




"sf 

r>- 

00 

o 

o 

\o 



00 

CO 

XT 

rH 

xr 

w 

to 

to 

iH 

o 

cn 

rx 



tH 

rH 

rH 



lO 


O 

f-H 

CO 

00 


• 


CO 

CX) 

to 

o 


rvi 


to 

,-H 

o 

Cl 


f 


rH 

rH 

rH 






xr 


(x 


to 


r--. 

CO 


O 


I 


to 

r-^ 


CJ) 




rH 

rH 




in 


o 

r'O 

CO 

rxj 


« 


to 

CO 

T— 

(XJ 




to 

r-^ 

o 

cn 


« 


f-H 

rH 

rH 





CjO 

vO 

(XJ 

o 


o 


to 

CO 

VO 

VO 


r— ( 


to 

rH 

O 

cn 


* 


rH 

rH 

rH 



lO 


\o 

to 

rx 


■■ 

" 


LO 

o 

cn 





to 

CX4 

o 



1 


r-t 

rH 

Bil 


Bl 



Ct> 

OO 

O 


bI 

to 


LO 

rH 

r-4 

to 


r—i 


to 

CXI 

rH 



* 


rH 

1 — 1 

rH 


■11 

U1 


CO 

CO 

OO 





1'- 

VO 

CO 





to 

Csl 

rH 



1 


rH 

rH 

■911 

PSI 




cn 

CX3 


CO 


C-J 


(^4 

CXI 


cn 


CNJ 


'<T 

to 


fH 


‘ 


rH 

rH 

m 

SHI 


to 


to 

rH 

M 



* 


00 

rH 




C'J 


XT 

XT 

ISB 



1 


rH 

rH 

Bl 



in 


to 

CO 

rH 

(x| 



vO 

XT 

rr 

to 




i/) 

to 

in 

to 


' 


rH 

rH 


•H 


m 


-rf 

to 

cn 

o 


' 


C'JI 

o 

CO 

c^ 


CH 


r-^ 

rH 

o 

o 


1 


Csl 

oi 

fvj 

f J 




r— t 

"f 

r--. 

c 1 


CD 


o 

rH 

o 

o 


to 


ai 

Oi 

cn 

cn 


1 


r-i 

('I 

C'4 

(XJ 









5 


Q 

Q 

ro 

O 


u 






cC 

to 

n: 

O 

p 

P 

O 









tn 






< 

ru 



«>/» 

nV*> 

ti > 

■'C 



o 

o 

O 

o 


f- 


•M-r- 

to 

1 

rH 





Scan No. 


in 


in 

P-3 

UJ 

a 

03 

Q 




00 

04 

On. 

VO 


o 

1 — t 


iH 

fH 

04 


JO 

c> 

CTt 

cn 

on 

cn 



rvi 

rvj 

04 

04 

04 




O 

o 

O 

On 





rH 

04 

rH 




O'. 

o 

cn 

cn 




r 4 

04 

04 

04 




f'- 

o 

cn 

tn 


cn 


o 

tn 

tn 

XT 


fsj 



r-r 

O'- 

r-N. 




H 

rH 

tH 

tH 


m 


o 

tr 

o 

XT 


• 


in 

O 

OO 

tn 


r-j 


un 

tn 

XT 

•4- 


r-Ji 


r~J 

rH 

H 

pH 

PIEI 



fNJ 

vO 

o 



O 



OO 

on 



rsi 



no 

on 





»H 

H 

H 


BE 

t-O 


^s| 

CO 

rH 

o 

miB 

• 



tn 

TT 

fN- 


r- 



KJ 

CH 

tH 


r—H 


r— H 

H 

fH 

rH 

hBb 



fVl 

00 

CO 

VO 


LO 


ro 

o 

OO 

cn 


rM 


-rr 

m 

rH 

o 




(H 

iH 

, — ( 

iH 


vn 


rH 

tn 

04 

tn 


• 


rH 

00 

tn 

XT 


fSl 



04 

rH 

o 


r— I 


r-H 

H 

tH 

rH 

IjHBBBI 



o 

On. 

CO 

o 


O 


o 

On. 

04 

cn 




xS 

04 

iH 

cn 




rH 

H 

tH 


eHeI 

tn 


pH 

O 

CO 

tn 

|B| 


r*4 

On 

OO 

On. 




*cr 

04 

O 

cn 




rH 

H 

H 







C ' 

cn 


H 

lO 


oo 

o 

Tn. 




XT 

CnI 

o 





f-H 

H 

H 


Bil 

LO 


rH 

cn 

m 

Xj- 



pH 

XT 

On 

04 


rvi 


'cr 

ONI 

o 

cn 



iH 

r-1 

fH 




un 

VO 

rH 

On. 

04 

tn 

O 

o 

fvj 

tn 

O'. 

04 

VO 

lO 

rr 

04 

o 

cn 

On 


Oi 

rH 

H 

pH 



urj 


rH 

vO 

^n 

CO 


• 


O 

to 

vO 

rH 


fNl 


XT 

04 

o 

cn 


1 


rH 

iH 

pH 





a> 

•O’ 

O 

vO 


in 


Cr> 

to 

tn 

04 


1 


tn 

04 

o 

cn 




rH 

fH 

r-H 



U*J 


fH 

O 

XT 

04 

■j^l 

« 


OO 

to 

in 

XT 




rO 

04 

o 



f 


rH 

H 

rH 


BUI 



O 

On. 

04 

CO 


o 


CX) 

XT 

On. 

XT 


rH 


tn 

04 

O 

cn 




rH 

rH 

r~( 


HH 

in 


o 

On. 

On 

rH 


rvi 



tn 

on 

cn 




04 

O 

cn 


rH 

1 


iH 

rH 

.H 



LO 


CTl 

^n 

VO 

CO 




00 

on 

m 


fH 


i~o 

04 

rH 

o 


' 


rH 

fH 

iH 

rH 


in 


XT 

00 

tn 



* 


O 

o 

CO 



r>» 


Tf 

tn 

iH 



r-H 

t 


«H 

»H 

tH 





O 

on 

on 

Rl 


o 


H 

On. 

in 



OJ 



tn 

04 

El 


1 


H 

rH 

pH 



un 


O'J 

On. 

O 



CNJ 


o 

tn 

xt 

XT 

CO 

tn 



o«) 


pH 

fH 

tH 

El 


in 


tn 

on 

o 




cn 

OO 

tn 



r^J 


in 

tn 

tn 



1 


H 

r-» 

rH 



in 


vO 

04 

cn 



r-^ 


tn 

o«. 

rH 

El 



fH 

rH 

rH 

El 


CM 


fvJ 

04 

<Ni 

Bl 




xr 

tn 

cn 



o 


XT 

VO 

XT 

Km 


tn 


cr» 

cn 

cn 

El 


> 


04 

04 

04 





Q 

Q 

Q 

Q 

oc 




VH 


VH 

cJ 

oC 

P 

P 

P 

P 

s 

3c 

O') 

< 

> 

r> 

pr*' 

>* 

UJ 



«?<*» 


tfW* 

Pw 


o 

o 

o 

O 




xr 

K) 

04 

rH 



















































































Scan No . 


VO 


CO 

W 




o 

w 

p 


o 

to 

12900 

1 

in 

CO 

CO 

rsj 

in 

o~ 

CO 

OJ 

2882 

2S79 


to 



O 




• 



o 

P" 

00 


r-- 



CO 

CO 

oo 





CM 

PM 

O.J 




cn 

cr> 

03 

03 


to 


rO 

CM 

CM 

CM 


CNJ 


r-- 

o. 

P" 

O' 




rH 

rH 

rH 

r-^ 


LO 


r—l 


CM 

to 




rO 

o- 

in 

PM 


rvi 


in 

•cr 

TT 

-?r 


C'J 


rH 

r■r^ 

fH 

tH 




O 

R 


CO 


o 


lO 

■SI 


to 


rsi 


"<cr 

■SI 


PM 




rH 

19 


f-H 


Ul 



R 

04 

VO 

imin 

m 




rH 

-tr 






CM 

rH 


r-H 



R 

rH 

i-H 

IH■| 



03 

•*3- 

03 

O 


in 


o 

to 

in 

CO 


rH 


'cr 

04 

rH 

O 




1 — ( 

rH 

rH 

rH 


LO 


ot 

00 


rH 

nn 

• 


\o 

rH 


rH 


rvj 


to 

CM 


o 


t— t 


rH 

rH 


rH 

mui 




•<T 

in 

CM 


o 



O'. 

vO 

O' 


rH 



rH 

O 

cn 






rH 





Ol 

cn 

to 

o 


Ul 


o- 

VD 

to 

to 


* 


to 

rH 

o 

on 




pH 

rH 

rH 


■HU 



cr\ 

XT 

O' 

oo 

Rjllll 



o- 

in 

CM 

o 


to 


to 

rH 

O 

cn 




rH 

iH 

rH 


BH 



cr> 

to 

03 



LO 


lO 

in 

03 

03 


• 


to 

iH 

C73 

CO 


<Ni 


rH 

rH 





m 

O^ 


\o 

VO 

O' 



\0 

rr 

O 

CO 

to 

o 

in 

K) 

rH 

o 

CO 

O' 


OJ 

rH 

rH 

rH 



LO 


rj- 

to 

iH 

O' 




O 


o 

OO 


(NJ 


to 

rH 

o 

CO 


1 


rH 

rH 

rH 





rH 

rH 


rH 


in 


vO 

H" 

03 

OO 


I 


to 

rH 

cn 

CO 




rH 

f-H 




tn 


03 

in 

o 

P'' 




XT 

XT 

rH 

O 


r-- 


to 

iH 

O 

03 


1 


tH 

rH 

rH 





o 

"XT 

vO 

CO 


o 


•rjr 

VO 

CM 

1 — t 


rH 


to 

rH 

O 

on 


' 


rH 

rH 

rH 



Ul 


OJ 

■XT 

CM 

O' 


• 


in 

O' 

tn 

vO 


C^l 


to 

rH 

o 

cn 


rH 


rH 

rH 

rH 





in 

rj- 

o 

cn 


t-O 


in 

03 

cn 

O 


rH 


to 

rH 

o 

O 


> 


fH 

rH 

r~H 

rH 


tn 


rH 

t'- 

cn 

cn 


* 


CTfc 

to 

i/) 

P' 


r''* 


to 

Ol 

rH 

o 


rH 

I 


rH 

rH 

r-H 

rH 




00 

\o 

O' 

PM 


O 


CM 

rH 

CM 

O' 


TM 


xf 

to 

CM 

rH 


t 


H 

rH 

1 — 1 

rH 


in 


cn 

to 

VO 

CO 





rH 

to 

O 


r^i 


"tJ* 

xr 

to 

ff-j 


PH 

» 


rH 

rH 

rH 

rH 




yO 

cn 

■rf 

rH 


in 


vO 

m 

to 

PM 


r>4 


to 

in 

to 

to 


1 


rH 

rH 

r— H 

rH 


m 


to 

rH 

rH 

CM 


• 


C3 

CM 

03 

o 


!>• 


1 — ( 

rH 

o 

rH 


04 

I 


03 

CM 

PV 

0*4 




o. 

VO 

cn 

VO 


O 


OV 

1 rH 

o 

O 


to 


00 

i CD 

cn 

03 


1 


OJ 

OJ 

OJ 

OJ 




Q 

a 

O 

Q i 


U 


1 — 1 



i—4 

cc 

UJ 

oc 

s 

o 

P 

a 








r- 

u\ 

< 

„ 

' 

1 


5 































no o o 


APPENDIX B 


C PROGRAM TO COMPUTE LINE INTEGRAL OF A GIVEN FUNCTION OF 

C TWO VARIABLES 

C DEFINITIONS : 

C A=LOWER LIMIT 

C B=UPPER LIMIT 

C D=MUMBER OF DIVISIONS 

C M=SLOPE OF THE GIVEN LINE 

C C=Y INTERCEPT 

C M=NUMBER OF READINGS FOR EACH VALUE OF THE SLOPE M 

REAL M 

0PEN(UMIT=1,FILE=’SUP3.DAT’ ) 

READ( 1 ,*)A,B,D,M,C 

OPEN (UNIT=H,FILE=’ SUP30.DAT' ) 

WRITE(5, 10) 

WRITE(2, 10) 

10 FORMAT ( 19X, ’THETA’ , 15X, ’S’ , 13X, ’ PROJ . DATA ’ /// ) 

PI=4-*ATAN{1. ) 

C=C-. 1 

H=(B-A)/D 

T=ATAN(M) 

T=T-PI/18. 


C .Compute for different values of M. 

DO 500 1=1,18 

E5- T=T+PI/18. 

R=T*180. /PI 
XA=ABS(R-a9.99) 

C BRANCH FOR THE CONDITION WHEN RAYS- ARE PARALLEL TO THE Y-AXIS 

IF(XA.LT.1. )CALL SUB1 (R,T, A,B,D,M,C) 

C=-1. 1/C0S(T) 

M=TAN(T) 

S=-1.1 

DO 500 N=1,21 
S=S+0. 1 

THE • IF-THEN-ELSE' STRUCTURE HELPS IN OVERCOMING THE ERROR 
IN THE TWENTY-FIRST ITERATION. 
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IF(M.LE.20)THEN 
C=C+0. 1/C0S(T> 

ELSE 

C=C+. '10001/C0S{T) 

END IF 

a=b-d*h 

S1=F(A,I'1,C)+F(B,M,C) 


C Compute the number of times the do loop is to be executed. 

C 

K=NIMT(D/H.-1. ) 

C 

C COMPUTE THE SUM OF ODD TERMS 


C 

so=o. 

A=A+H 

DO 100 KI=1,K 
P=^F(A.M,C) 

SO=SO+P 

A=A+S.*H 

100 CONTINUE 

C INITIALIZE A FOR COMPUTING THE SUM OF EVEN TERMS 

A=B-(D~2. )*H 


C 

C 

C 


200 

C 

C 

C 


500 

400 

C 

C . 
C 


COMPUTE THE SUM OF EVEN TERMS 


BE=0. 

DO 200 KJ=1,K 
Q=F{A.M,C> 

SE=SE+Q 

A=A+2.*H 

CONTINUE 

APPLY SIMPSON’S ONE-THIRD RULE 

DI=(H/3. )*(S1+i4.*S0)+C2.*SE) )*(SQRT ( 1 .+M**2) ) 

R=T«iaO./Pl 

WRITE(5.400)R,S,DI 

WRITE(2,400)R,S,DI 

CONTINUE 

STOP 

FORMAT ( 13 X,F! 5 .S, 3 X,F 12 . 2 , 7 X,F 13 .a) 

..End of main program. 

END 


C 

C 

c 
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o n o 


THIS FUNCTION DEFINES THE VALUE OF F(X,Y) ALONG ALL LINES 
OTHER THAN THE Y-AXIS AND LINES PARALLEL TO IT. 

c 

FUNCTION F(X,M,C) 

REAL H 
Y=mX■^C 
SQ=y*i{g4- X*^*2 
IFCSQ.LT. 1. )THEM 
F=1. 

ELSE 

F=^0. 

EMDIF 

RETURN 

E^4D 

C 

C 

C THIS FUNCTION DEFINES THE VALUE OF F<X,Y) ALONG THE Y-AXIS 

C AND LINES PARALLEL TO IT. 

C 

C ^ -Jf 

C 

FUNCTION G(li,C,Y) 

REAL M 

SQ1-Y^^t2+C-M-Jf^2 
IF(SQ1.GE.1. )THEN 
G=0. 

ELSE 

G=1. 

EMDIF 

RETURN 

END 

SUBROUTINE SUB1 (R,T, A,B,D.N,C) 

REAL M 
C 

PI=4-^f-ATAN( 1 . ) 

C 

C THE FUNCTION DEFINITION ENSURES THAT , IF DESIRED, INTEGRALS OUTSIDE A 

C UNIT CIRCLE ARE ZERO. 

C 

C WE HAVE TO RESORT TO A SUBROUTINE , BECAUSE FOR THE Y-AXIS, 

C THE SLOPE 'M* IS INFINITY. THE ESSENTIAL DIFFERENCE BETWEEN 

C THE SUBROUTINE AND THE MAIN PROGRAM IS THAT, FOR COMPUTING 

C THE LINE INTEGRAL, ONLY THE LIMITS ON Y NEED TO BE KNOWN 

A=-1.0 
B=1 .0 
H=(B-A)/D 
C=-1.1 


65 



DO 501 N=1,E1 
JF(M.LE.EO)THEN 
C=C+. I 
ELSE 

C=C+. 10001 

ENDIF 

A=B-DflH 

S1=G(M,C,A)+G(M,C.B) 

C 


C Compute the number of times the do loop is to be execute 

C 

K=NIMT(D/E.-1 , ) 

C 

C COMPUTE THE SUM OF ODD TERMS 


C 

so=o . 

A=A+H 

DO 101 KI=1,K 
P=G(M,C, A) 

SO=SO+P 
A^A+E . ■H'H 

101 CONTINUE 

A=B- ( D-S . ) -K-H 
C 

C COMPUTE THE SUM OF EVEN TERMS 

C 

SE=0. 

DO SOI KJ=1,K 
Q=G ( M , C , A ) 

SE=SE+Q 
A=A+E . 

SOI CONTINUE 

C 

C APPLY SIMPSON’S ONE-THIRD RULE 

C 

DI=(H/3. )*(S1+(4.*S0)+(S.*SE) ) 

R=90» 

WRITE(5,400>R,C,DI 
WRITE(E, 1-00)R,C,DI 
501 CONTINUE 

C 

C UPDATE THE VALUE OF THE ANGLE TO R=100. DEGREES. 

C 

T=T+PI/18. 

400 F0RMAT(13X,F15.8,3X,F1S.S,7X,F13.8) 

RETURN 

END 
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APPENDIX C 


C !HI -it H-H #-}f ■}! •« •«•■«•■« ********' 

C PROGRAM FOR FAN-BEAM DATA. 


C P(41 )=MATRIX OF THE PROJECTED DATA 

C FAN(756,3>=MATR1X OF FAN-BEAM DATA CONVERTED TO PARALLEL 

C Q{41)=MATRIX OF FILTER VALUES 

C Q1(41 ,41 )=MATRIX OF ’Q’ VALUES PROJECTED DIAGONALLY 

C PIC1 (41 )=MATRIX OF CONVOLVED VALUES 

C PIC(41 ,41 )=MATRIX OF VALUES FROM PREVIOUS SCAN 

C FPIC(41,41>=FINAL PICTURE 

C IPIC(41 ,41 )=INTEGER VALUES OF FPIC FOR SCREEN DISPLAY 


DIMENSION FANC 756,3) ,P (41 ) , (3 ( 4 1 ) , Q1 ( 41 , 41 ) , PIC ( 41 , 41 ) , PIC 1 ( 
DIMENSION IPIC(41,41) ,FPIC ( 4 1 , 4 1) , TEMP ( 3 ) ,PRQJ (378,2) ,DEN(2 


C RAM.FIL CONTAINS THE '(3' VALUES 

C PIC. DAT STORES THE VIEW AT A GIVEN ANGLE 

C FPIC.DAT STORES THE IMAGE AFTER EVERY ITERATION 

C DIV.DAT CONTAINS THE NUMBER OF PIXELS IN A ROW 

C P1D.T10 CONTAINS PROJECTED DATA 

C AIR. REF CONTAINS THE DATA FOR PLEXIGLASS CORRECTION 

C***************-^*}* **'«-'}t-*'H'*-»-«-**-«'**-«-'«-**#**#-«--»--M-*'Jf-**-«-**-»*******'**-**-*t-****' 


OPEN(UNIT=EO,FILE=’RAM.FIL’ ) 
0PEN(UNIT=E2,FILE=''DIV1D.T10’ ) 
0PEN(UNIT=E6,FILE='P1D.T10’ ) 
OPEN ( UNIT=E3 , FILE= ’ DI V . DAT ’ ) 
OPEN ( UNIT=28 , FILE= ' AIR . REF ’ ) 


PI=4 .*ATAN( 1 . ) 

C READ INPUT PARAMETERS 

READ ( 23 , ■»• ) MDI V , MVIEW , D , RAD 
SGMAX=D*SIN(25-*PI/180. ) 

DO 75 1-1,378 

READ (26,^0 BETA , SGMA , DATA 
FAN (1,1) -DETA+SGMA 

FAN ( 1 , E ) =D-3!-SIM ( SGMAH-P 1/180.)/ SGMAX 
FAN (1,3) =DATA 
75 CONTINUE 
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C INTERPOLATION 

1=1 

DO 80 J=1 , 18 
XX=-1 .0 
PROJ (I, i )=XX 
PROJ (I,2)=FAN( 1,3) 

DO 80 K=1,E1 
1 = 1+1 

IF(I.GT.37S)G0T0 65 

XX=XX+.1 

XX1=FAM(I-1,2) 

XXE=FAM(I,2) 

YY1=FAN{I-1 ,3) 

YY2=FAN(I,3) 

PROJ (1, 1 )=XX 

PRaj(I,2)=YY1-( (XX 1 -XX)*(YY 1 -YY 2 )/(XX 1 -XX 2 ) ) 

80 CONTINUE 

C READ THE DATA FOR PLEXIGLASS CORRECTION INTO AN ARRAY 

65 READ(28,51 ) ( (DEN<M,N) ,M=1 ,21 ) ,N=1 ,21 ) 

C 

C INITIALIZE ALL PIXELS TO ZERO 

C 

DO 500 I=1,NDIV 
DO 500 J=1,NDIV 
FPIC(I,J)=0. 

PICd, J)=0, 

IPIC( l,J )=0 
500 CONTINUE 

C READ IN THE VALUES OF THE Q FUNCTION, FOR CONVOLUTION 

READ(EO, 11 ) {Q(N) ,N=1 ,NDIV) 

C DISPLAY FORMAT MESSAGE 

WRITE{5, 10) 

WRITE (22, 10) 

C 

C SPREAD ELEMENTS OF Q ALONG DIAGONALS OF Q1 

C 

DO 101 N1=1,NDIV 
DO 101 i11 = 1,NDIV 

NMD=ABS(N1-I'11 ) + 1 
01(N1,M1 )=Q(NMD) 

101 CONTINUE 

C 

C START VIFWIMG AT DIFFERENT ANGLES 

C 

X=0. 

Kl<=0 
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non W non 


DO 700 1=0, (NVIEW-1 ) 

X=X-H . 

13=0 

DO ai II=UMDIV 
13=13+1 
|{K=KK+ 1 

Pdl )=-L0G(PR0J (13,2) ) 
81 CONTINUE 


CARRY OUT CONVOLUTION 

DO 200 II=T,NDIV 
PIC 1(II)=0. 

DO 200 J !=1 ,NDIV 

PIC1(II>=PIC1(II) + ((31(II,J1)^P(J1) ) 

CONTINUE 

T=REAL { I ) +PI/REAL( NVIEW ) 


ADJUST FOR OVERLAP FOR DIFFERENT VIEWS 


NDV=NDIV/2+ 1 
NDV1=NDV-1 

DO 699 ri1=-NDV1 ,NDV1 

DO 699 N1=-NDV1 ,NDV1 
N=M1+NDV 
N=N1+MDV 
Rri=REAL(M) 

RN=REAL(N) 

DGQ=ABS ( S(3RT ( ( RM-NDV ) *^<2+ ( RN-NDV ) «*2 ) ) 

IF(RM.EQ.NDV)THEN 

PHI=PI/2. 

ELSE 

PHI=ATAN ( ( RN-NDV ) / ( RM-NDV ) ) 

END IF 

R= ( DS(5^« COS ( T-PHI ) ) +NDV 
L=INT(R) 

1F(L.LT. 1.0R.L.GT.20)G0T0 699 

S=(R-REAL(L) )*(PIC1 (L+1 )-PIC1 (L) )+PIC1 (L) 
Z=(E-0*PI*<REAL(NVIEW) )*REAL(NDIV-1 ) ) / ( REAL ( NDIV**2 ) 
PIC(M,N)=FPIC(M,N) 

FPIC (M,N) = ( (FPIC (M,N)-»(X-1 . ) ) + ( ( S*Z-DEN ( M , N ) ) /RAD) ) /X 
C EQUATE THE PIXEL VALUES OUTSIDE CIRCLE TO ZERO 
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MSQ=-- ( M-NDV ) ( N-NDV ) **2 

IF C MSQ . GT . NDV 2 } THEN 
FP1C(M,N>=0.00 
ELSE 
ENDIF 

IPIC(i'1,N)=ABSCNINT(FPIC(M,N)#100. ) ) 

699 CONTINUE 

C THIS OPTION ALLOWS US TO TERMINATE THE PROCEEDINGS IF THE ERROR OF 

C RECONSTRUCTION REACHES A PRESET LOWER LIMIT 

SDIS=0. 

SPIC=0. 

DO 695 M=1,NDIV,3 
DO 695 N=1,NDIV,3 

DIS=ABS( (FPIC(M,N)-PIC<M,N) ) ) 

SDIB=SDIS+DIS 
SPIC=SPIC+ABS(FPIC{M,N) ) 

695 CONTINUE • , 

IF<SDIS-LT. ( .005*SPIC) >GOTO 694 

DEG=iaO.«-T/PI 

WRITE (5, 50) ( (IPIC{M,N) ,N=1,NDIV) ,M=1,NDIV) 

WR1TE(5,111)DEG 

WRnE(21,1-11)DEG 

50 F0RMAT(21(6X, 21(13)/)/) 

51 F0RMAT(21(8X,F14.8/)/) 

15 FORMAT ( ax, A72///) 

100 FORMAT ( 50X, F I3. 8) 

111 FORMAT (fl5X,F10. 5, 2X, 'DEGREES’ ) 

700 CONTINUE 

694 WRnE(22,50) ( (IPIC(M,N) ,N=1 ,NDIV) ,M=1 ,ND1V) 

C CALCULATE THE MAXIMUM AMD MINIMUM PIXEL (LITF) VALUES 

RM3N=0. 

RMAX=--0. 

DO 800 I=1,NDIV 

DO 800 J=1,NUIV 

IF(FPIC{I, J ) .GT,RMAX)RMAX=^FPIC(1, J ) 

TF(FPIC<1, J ) .LT.RMIN)RMIN=FPIC (I, J ) 

800 CONTINUE 

WRHE(22,13)RtlAX,RMIN 

C CALCULATE THE AVERAGE PIXEL (LITF) VALUE 


SUM=0 . 
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mm ^ 


900 


10 

11 


DTR--=0 „ 

DO 900 I=3,NDlV-4 

DO 900 J =-3,11019-4 
SUfi^^SUMi-FPICd^J ) 

DTR=-DTR+1 . 

CONTINUE 

AVG=SUi1/DTR 

WRITE(H2, 12)NDIV,MDV 

KfR ITE ( cE ,5 1 ) ( FP IC ( 1 ! , N ) , N= 1 , MDI V ) 

WRITE (EE, 14)AVG 
WRITE(5,*)AVG,DTR 

WRITE (EE, 51 ) ( (FPIC(i1,N) ,M=1,NDIV) ,M=1 ,NDIV) 

FORMAT (EX d HORIZONTAL CENTERLINE ROW IE/ COLUMN 1 TO’,. 
FORMAT (EX, ’MAXIMUM LITF=’ , EX , F14 . 8 , EX , ’MINIMUM LITF=' ,F14.S 
FORMAT (EX, ’ AVERAGE= ’ , F14 . 8// ) 

STOF^ 

FORMAT (EX, ’POSITIVE IMAGE’//, EX, ’ FPIC . DAT ’ // ,2X , 'PAM.FIL’ ) 
F0RMAT(3(9X,E15.7) ) 

END 





APPENDIX 


D 


-CTO OEMERATE A SERIES RANDOM NUMBERS SUPERIMPOSE THEM 

DATA . > 

PROGRAM MOV? 

VAR 

R,SGRT i :REALj 
C , N , SRRT2 , MPOT : li'l ITGER ? 

IP, OP; TEXT; 

BEGIN 

ASS.1 GM ( IP , ’ C ; T !ODD „ DAT ’ ) ; 

RESET < IP ) ? 

ASSlGhKOP, 'C:rOl 1D.T10' ) ? 

REWRITE (OP) ; 

C s =• ! ,- 

WHILE V, <~-2i DO 
BEGIN 

READLM(1P,R,N) ; 

SQiU1 - =^SQRT(N) j 

sern 2 : ==round ( sort 1 ) ? 

NPOI ; =2*nAIMD0!1 ( SORTS ) -B(5RTE ; 

N!:^Mi-( l-HNPOI) 

WRITELN(OP,R5 10=3, N: 10) ? 

C;=C+-i; 

END: -rHHII..E> 

CLOSE COP) ? 

END. -fRANDNUrO 


TP 


0 ^! 



APPENDIX E 


RECONSTRUCTED INAGE FDR PINE. 


rOSITJVE I!1A(^E 
FPIC . DAT 
RAH.FIL 

0 0 0 0 0 0 000050000000000 

0 0 0 0 0 0 5 6 7 7 7 7 7 7 7 0 0 0 0 0 0 

0 0 0 0 4678688388777 0 0 0 0 

0 0 0 57 8 808880 8 88777000 

00478888888888887770 0 

0 0 6 8 5 8888888688887700 

0 3 7 8 86888999988888770 

0 6 838888999999888877 0 

078 8 88899999999888870 

0 7 8 886899999999888870 

5 7 8 688S9999999SS88875 

0 7 808899999 9 998888870 

0 76888999999998888870 

0 7 7 8 6 6899999988888860 

077088889999888888750 
0 07 7 0688888888888860 0 

0 0 7776888888883887400 

0 0 0 7 7 7 8 88388888875000 

0 0 0 0 7 7 7 8 8 8 8 8 8 8 7 6 4 0 0 0 0 

0 00 0 0 0 7 77777765000000 

0 00 0 0 0 000050000000000 


!1AXf!il'li LIir-=^ 0.06874133 lilNinUfl LITF= -0.05289834 

HUR 17-ONTAi.. CF.Ml ERL] ME 
R 01-121 

COLUMN ! TOli 

0,0500??. 77 
0„ 0704 9052 
0.077055 15 
0 . 07877i664 
0,06 131342 
0„()a379.398 
0.08466700 
0.037066 U3 
0,06687436 
0,06755? I I 
0,06603/49 
0.06755311 
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AVERAGE---^ 


0.0A687136 
0.03703610 
0.08463730 
0 „ 06379398 
0,08131342 
0.07876661 
0.07905513 
0.07049052 
0.05009277 


0.08329582 

0,00000000 
0 . 00000000 
0 . OOOOOOOO 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
O. OOOOOOOO 
0.00000000 
0.00000000 
0,05009277 
0.00000000 
0 . OOOOOOOO 
0 . OOOOOOOO 
0 . OOOOOOOO 
0.00000000 
O. OOOOOOOO 
0 . OOOOOOOO 
0.00000000 
O , OOOOOOOO 
0,00000000 


O „ OOOOOOOO 
0.00000000 
O , OOOOOOOO 
0 , OOOOOOOO 
0 . OOOOOOOO 
0 . OOOOOOOO 
0.04068446 
O .<)59ej3553 
0.06622 172 
0 . 06966973 
0,07049052 
0.06045242 
O . 06726465 
0,067000 12 
06966036 
o.oOOOOOOO 



0^00000000 
0.00000000 
O . 00000000 
0 . 00000000 
0 . 00000000 


0 00000000 
0 . 00000000 
0.00000000 
0.00000000 
0.0^300aZ5 
0.06 i679^6 
0.07 1E9 160 
0.076 I3724 
0.07 BE 1 389 
0. 07920318 
0.,()79055ir> 
0.078-14675 
0.07728843 
0.07486982 
0.06928518 
0.06817999 
0.07150093 
()., 00000000 
0,00000000 
0 , 00000000 
0.00000000 


0.00000000 
{). 00000000 
O . 00000000 
0.04754804 
0.06619010 
0.07559633 
0 . 07900287 
0.07B86074 
0,07843288 
0.07870436 
0.07876664 
0.07904380 
0.07807133 
0.07781857 
0-07748792 
0.07406854 
0.06/55796 
0.0/019676 
0 . 00000000 
0 , 00000000 
0.00000000 
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Date Slip 

This book is to be returned on the 
date last stamped. 
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